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Abstract — We examined the effects of inhomogeneity on the 
dynamics and structural properties using Boolean networks. Two 
different power-law rank outdegree distributions were embedded 
to determine the role of hubs. The degree of randomness and 
coherence of the binary sequence in the networks were measured 
by entropy and mutual information, depending on the number 
of outdegrees and types of Boolean functions for the hub. With 
a large number of outdegrees, the path length from the hub 
reduces as well as the effects of Boolean function on the hub 
are more prominent. These results indicate that the hubs play 
important roles in networks' dynamics and structural properties. 
By comparing the effect of the skewness of the two different 
power-law rank distributions, we found that networks with 
more uniform distribution exhibit shorter average path length 
and higher event probability of coherence but lower degree of 
coherence. Networks with more skewed rank distribution have 
complementary properties. These results indicate that highly 
connected hubs provide an effective route for propagating their 
signals to the entire network. 

I. Introduction 

Recent studies of complex networks [1] have shown inho- 
mogeneous connectivities having a small number of highly 
connected nodes, frequently called as hubs, along with many 
poorly connected nodes. The inhomogeneity in connectivity 
has large effects on the property and/or function of the hubs; 
for example robustness in metabolic, genetic regulatory and 
neural networks [2] ; performance of artificial neural networks 
. [3]; and dynamics of coupled oscillators [4], [5]. 

The Boolean network [6] is one of the discrete dynam- 
ical models for the transcriptional regulatory network and 
exhibits binary sequences of the state variables that represent 
expression pattern of the network [6], [7]. Since the state 
variables in the network are sensitive to inputs from other 
nodes via directed edges, and affect other nodes, the quality 
of communication is characterized as the size of mutual 
information [8]. This mutual information indicates the degree 
of coherence, synchronization, amount of information content 
in the state variables, or potential for computational capability 
of the network [9]. 

In this study, we show the role of hubs for the emergence of 
coherence in a Boolean network. Since we embedded power- 
law rank outdegree distributions in the Boolean networks with 
an input connectivity of Ki n = 2, the model networks have 
some hubs that integrate many outdegree connections. Because 
the hubs synchronously transmit their state to the downstream 
nodes, these nodes are simultaneously affected by single or 
multiple hubs . The structural condition seems to automatically 
provide global coherence in the state variables; however, the 
structural aspects give only the possible effects of the hubs. 



In fact, we need to consider a type of Boolean function at 
the hubs and path length from the hubs. We show both the 
effects of Boolean functions and path length on both the event 
probability and size of entropy and mutual information. 



II. Model 

Dynamics of the Boolean networks [6], [7] are determined 



by 



Xi(t + l)=Bi[X(t)] (i = 1,2, ...,7V), 



(1) 



where Xi(t) is the binary state, or 1, of node i at time 
t\ Bi(-) is the Boolean function [see Tables I— III] used to 
simultaneously update the state of node i\ and X(t) is a binary 
vector that gives the states of the N nodes in the network [See 
Fig. 2(a)]. After assigning initial states X(0) to the nodes, 
the successive states of the nodes are updated by input states 
and their Boolean function. The dynamical behavior of these 
networks is represented by the time series of the binary states. 
The time course follows a transient phase from an initial state 
until a periodic pattern, known as an attractor, is established. 



TABLE I 

AND-OR type Boolean functions with indegree Ka* 



■■ 2. Since 



OUTPUT PROBABILITY OF ONE OR ZERO EQUALS \ OR | WITH RANDOM 



BINARY INPUTS THE FUNCTIONS HAVE THE SAME OUTPUT ENTROPY, 

(|) =#(!)= 0.811. 



Input 


Output 




1 

1 
1 1 


1 1 1 1 
10 10 11 

10 110 1 

1 1 1 1 



TABLE II 

XOR type Boolean functions with indegree Kn, 



■- 2. Since 



OUTPUT PROBABILITY OF ONE OR ZERO EQUALS ± WITH RANDOM BINARY 



INPUTS THE FUNCTIONS HAVE THE SAME OUTPUT ENTROPY, 



i. 



Input 


Output 




1 

1 
1 1 


1 1 1 
110 1 
10 10 10 
110 10 



III. Numerical Condition 

We randomly constructed 10 4 Boolean networks in each 
power-law rank distribution [see Fig. 1] with a fixed network 



TABLE III 

CONSTANT type Boolean functions with indegree K in = 2. 
Since output probability of one or zero equals 1 or with 
random binary inputs the functions have the same output 

ENTROPY, H(0) = H(l) = 0. 



Input 


Output 





1 


1 


1 


1 


1 


1 1 


1 



size. 2 xlO 3 initial states were applied to each network. 
Sixteen different Boolean functions [see Tables I— III] were 
used with equal probabilities. Note that all the generated 
networks use the same amount of resources since the size of 
the network is fixed-256 nodes, 512 directed edges, and 16 
sets (16 2 = 256) for all the Boolean functions [see Table VI in 
Section VI]. We measured entropy (randomness) and mutual 
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Fig. 1. Power-law rank connectivity distributions in the model. Power 
exponent, 7 for type I and type II is about 0.8 and 0.5, respectively, where 
K ou t(rank) ~ rank -1 . We performed only a single network size of 256 
(= N) in this paper. 



information (coherence) of the state variables to characterize 
the dynamics of the Boolean networks [7], [8], [9] [see Fig. 
2]. 



IV. Results 



A. Dynamics 



In total, we obtained 137254 and 459240 attractors from 
type I and type II distributions, respectively. The size of 
entropy and mutual information were measured from the 
attractors. Major statistics of dynamical properties are shown 
in Table IV. By comparing the effect of the skewness of 
the two different power-law rank distributions, we found that 
on average networks with more uniform distribution (type II) 
exhibit higher event probability of coherence, but lower degree 
of coherence. 



(a) 




(b) 













Fig. 2. (a) Flow of state variables from upstream to downstream. Input 
sequence Inn for node i is the same as the output sequence of an upstream 
node Outji, and the output sequence Outi for node i is the same as 
the input sequence of a downstream node Ink- When node i has multiple 
output connections, it has the same binary sequence because state variables 
in networks are subject to Eq. (1). (b) Typical local structure around a hub. 
Squares and directed edges (arrows) correspond to the nodes and connections 
of binary sequence pathways, respectively. H in the squares corresponds to 
the Boolean function. Each node has one of the Boolean functions as shown 
in Tables I— III [see also Section VI]. 



TABLE IV 

Dependence of dynamical properties on outdegree 
distributions: magnitude relations are indicated. 



Properties 


Type I 




Type II 


No. of attractors 


from 10 4 nets 


137254 


< 


459240 


Proportion of 


positive entropy 


0.788 


< 


0.848 


Proportion of 


positive mutual information 


0.661 


< 


0.698 




Median 


34.8 


> 


31.6 


Entropy / bit 


Average 


46.1 


> 


38.4 




3rd quartile 


83.3 


> 


63.2 


Mutual 


Median 


0.827 


> 


0.387 


information / bit 


Average 


2.48 


> 


0.996 




3rd quartile 


3.73 


> 


1.36 



B. Path length 

To get stractual property of propagating route of state 
variables, we measured two properties: 

1) Path length which the average number of the directed 
edges in the shortest path from a node to all reachable 
nodes. 

2) Average path length which the average number of the 
path lengths for all the nodes. 

Figure 3 shows the differences in the both path lengths. Nodes 
with higher outdegrees have shorter path lengths. 

Theoretical relationship between the number of outdegrees, 
X of a starting node and path length, i can be written as [For 
details, see Section VI.] 



/255 

1 = Log Kin ( — 



1 



(2) 



We obtained good relationship between numeical results and 
Eq. (2). 
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Fig. 3. (a) Relationships between the rank of the hubs [see Fig. 1] and the 
path length, (b) Dependence of path length on the number of outdegree from 
a starting node. The dotted lines show the relationship of Eq. (2) [see Section 
VI]. Both path lengths are obtained from 10 4 generated networks with TV = 
256. 



C. Rank Dependent Dynamics 

Together with the network structural condition, we show 
the dependence of rank distributions on entropy and mutual 
information in Figs. 4 and 5. We classified different 16 
Boolean functions into three types (AND-OR, XOR, and 
CONSTANT types) based on the input-output relationships 
[10], [11]. The dependence of the Boolean function on entropy 
and mutual information is prominent in higher ranked hubs on 
both the outdegree distribution styles. The Boolean functions 
with larger entropy tend to have larger mutual information. 
These results have weak dependence but clear tendencies, 
and suggest that the collective (global) coherence in the state 
variables of the networks is subject to the style of upstream 
(local) conditions, including the number of outdegrees from a 
hub and assignment of the Boolean functions. 

V. Conclusions 

We have summarized the results in Table V. A critical 
control parameter for synchronizability of coupled oscillators 
in complex networks is coupling strength or path length among 
the oscillators [4], [5]. Since our Boolean networks consist of 
the same amount of resources, some nodes integrate many 
edges, while others have less. As shown in Fig. 3, with a 
larger number of outdegree from hubs, the path length from the 
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Fig. 4. Dependence of proportion of positive entropy and size of entropy (a) 
and (b) on type I and type II distribution in Fig. 1. Three different symbols 
denote different types of Boolean functions of the hubs: square; CONSTANT 
type; triangle; AND-OR type; and circle; XOR type [see Tables I— III] . H 
denotes the Boolean functions of the hubs [see Fig. 2 and Table VI in Section 
VI]. 



hubs reduces. The shorter path length can result in a stronger 
coupling strength between the hubs and other nodes. Because 
of the differences between type I and type II distributions, the 
higher ranked hub in type I distribution have an advantage to 
transmit their signals [see Figs. 4(a) and 5(a)], and the hubs in 
type II distribution have longer path lengths and show smaller 
dependences [see Figs. 4(b) and 5(b)]. On the other hand, 
type II distributions exhibit higher event probability of entropy 
and mutual information and a large number of attractors [see 
Tables IV and V]. The resultant dynamical properties are 
due to the more decentralized topology that can lead to the 
presence of some mid-rage hubs. The hubs permit local and/or 
weak coherences to emerge in the networks. 

The higher indegree Ki n > 2 provides higher density of the 
edges as well as a narrower range of the path length [see Fig. 
3], and therefore the dependence of the outdegree from a hub 
will become less prominent. 
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Fig. 5. Dependence of proportion of positive mutual information and size 
of mutual information (a) and (b) on type I and type II distribution in Fig. 1. 
Three different symbols denote different types of Boolean functions on the 
hubs: square; CONSTANT type; triangle; AND-OR type; and circle; XOR 
type [see Tables I— III] . H denotes the Boolean functions of the hubs [see Fig. 
2 and Table VI in Section VI]. 



VI. Appendix 

Since the networks are randomly constructed with K in = 
2, all the nodes have two downstream nodes, on average, 
regardless of outdegree distribution styles. The total number 
of downstream nodes from a starting node T(i) can be written 

as 



T(i) = ^2 n = 2 l -1, 



(3) 



where i is the number of steps (path length) from the starting 
node. With increasing i, the total number of downsteam nodes 
increases. When the starting node has X outdegree, 



XT(i) =X(2 i - 1). 



(4) 



Since the size of network is limited to 256 in the paper, the 
total number will reach to the 255. 



X(2 l - 1) = 255, 



(5) 



where 255 indicates the exclusion of a starting node from the 
network. Finally, we can obtain the number of steps (path 



TABLE V 

Summary of results: Magnitude relations are indicated. 



Dynamical and structural properties 


Types 


Boolean function dependence 


XOR>AND-OR>CONSTANT 


Skewness of outdegree distribution 


I > II 


Average path length 


I > II 


Size of entropy and mutual information 


I > II 


Rank dependence 


I > II 


Event probability of positive 
entropy and mutual information 


I < II 



length), 



'255 

i = Log 2 ( — 



1 



(6) 



When X = 1, the i takes 8, meaning that 8 steps are required 
to reach other nodes in the same network on average [see Fig. 
3(b)]. For general form of Eq.(6), 

'(N-l)(K in -l) , 



i Log Ki , 



X 



(J) 



where Ki n is indegree of the Boolean network and N is the 
number of nodes in the network. 



TABLE VI 

Number of networks out of 10 4 networks (type I) are 
indicated. Three different types of Boolean functions [see 
Tables I-III] at the hubs are used with equal probability. "H" in 

THE TABLE DENOTES THE HUBS BOOLEAN FUNCTION [SEE FlG.2(B)]. 



Rank 


H=XOR 


H = AND-OR 


H=CONSTANT 


1 


3725 


4995 


1280 


2 


3736 


5042 


1222 


3 


3696 


5029 


1275 


4 


3710 


5020 


1270 


5 


3705 


5095 


1200 


6 


3725 


4996 


1279 


7 


3756 


5006 


1238 


8 


3768 


5012 


1220 


9 


3815 


4951 


1234 
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